Accepted in Astrophys. J. 

Preprint typeset using I4T]^X style emulateapj v. 03/07/07 



EXPLOSIVE NUCLEOSYNTHESIS IN THE NEUTRINO-DRIVEN ASPHERICAL SUPERNOVA EXPLOSION OF A 

NON-ROTATING 15M STAR WITH SOLAR METALLICITY 

Shin-ichiro Fujimoto 1 , Kei Kotake 2 , Masa-aki Hashimoto 3 , Masaomi Ono 3 , and Naofumi Ohnishi 4 

Accepted in Astrophys. J. 



o 

(N 

C 
3 



6 

CO 

o3 



> 
O 
(N 

O 



13 



ABSTRACT 

We investigate explosive nucleosynthesis in a non-rotating 15M Q star with solar metallicity that explodes 
by a neutrino-heating supernova (SN) mechanism aided by both standing accretion shock instability (SASI) 
and convection. To trigger explosions in our two-dimensional hydrodynamic simulations, we approximate 
the neutrino transport with a simple light-bulb scheme and systematically change the neutrino fluxes emitted 
from the protoneutron star. By a post-processing calculation, we evaluate abundances and masses of the SN 
ejecta for nuclei with the mass number < 70 employing a large nuclear reaction network. Aspherical abun- 
dance distributions, which are observed in nearby core-collapse SN remnants, are obtained for the non-rotating 
spherically-symmetric progenitor, due to the growth of low-mode SASI. Abundance pattern of the supernova 
ejecta is similar to that of the solar system for models whose masses ranges (0.4— 0.5)Mq of the ejecta from the 
inner region (< 10, 000km) of the precollapse core. For the models, the explosion energies and the 56 Ni masses 
are ~ 10 51 erg and (0.05-0.06)M Q , respectively; their estimated baryonic masses of the neutron star are com- 
parable to the ones observed in neutron-star binaries. These findings may have little uncertainty because most 
of the ejecta is composed by matter that is heated via the shock wave and has relatively definite abundances. 
The abundance ratios for Ne, Mg, Si and Fe observed in Cygnus loop are well reproduced with the SN ejecta 
from an inner region of the 15M progenitor. 

Subject headings: Nuclear reactions, nucleosynthesis, abundances — stars: supernovae: general — Hydro- 
daynamics — Methods: numerical 



1. INTRODUCTION 

The explosion mechanism of core-collapse supernovae 
(SNe) is still not clearly understood. Multi-dimensional ef- 
fects such as standing accretion shock instability (SASI) and 
convection are recognized to be most important for unveiling 
the explosion mechanism, in particular for a prog enitor heav- 
ier than about 1 IMp, in it s main sequence phase dKitaura et alj 
l2006tlB~uras et alj2006alfbh . Here SASI, becoming very popu- 
lar in current supernova researches, is a uni- and bipolar slosh- 
ing of the stalled supernova shock with pulsationa l stron g 
expansion and cont ra ction (see, e.g . lBlondin et al.1 8 2003): 
ISchecketaT] (120041): fOhnishi et alfoOOd): iFoglizzo et alj 
(l2007l):lBlondin & M ezzacappal (l2007h : llwakami et al J (I2008L 
2009); Nordhau s et al.l (12010) and references therein). Some 
of recent two-dimensional (2D) radiation-hydrodynamic sim- 
ulations show that the delayed neutrino-driven mechanism 
aided by SASI and convection does work to produce as- 
pherical explosio ns dMarek et aT1l2009l:lM~arek & Jankall2009l: 
ISuwa et aUhoIflh . 

Observationally, global anisotropics and mixing as well 
as smaller-scale clumping of the SN ejecta, a re co mmon 
featur es of S N remnants like in SN1987A dWang * 

2002), Cas A dHughes et "all 1200(3: IWillingale et aljj l 

G292 0+1.8 dPark et all 120071). and Cygnus loop dKjmural 
120091: lUchida et al.ll2009l) . Asymmetries commonly observed 
in the nebular emission-line profiles are considered as an 
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evidence that core -collap se SNe occur gene r ally aspheri- 
callv dMaeda et all 1200a [Modiaz et all 120081: iTanaka et alj 
2009b; TaubenbergeretdJll009). Evidences for asymme- 
try are also obtained from spectropolarimetric obs ervations of 
Type Ibc SN a t an early phase (~ days) (see, e.g..Tana kaet alj 
(2008, 2009a) and references therein). 

Thus far, nucleosynthesis studies of the SN ejecta have 
almost successfully reproduced the solar comp osition and 
abund a nces of radioactives o bserved in SN1987A dHashimotol 
1995b I Woosley & Weaver! 11995b iThielemann et all 119961: 
Rauscher et al.l l2002j). However those spherical models have 



some problems such as overproduction of neutron-rich Ni 
isotope s and underproducti ons of 44 Ti, 64 Zn and light p- 
nuclei (Rauscher et al. 2 0021) . 

Aspherical effects on the explos i ve nu c leosynthesis hav e 
been investigated bv lNagataki et alj d 1 9971) : iNagatakil d2000t) . 
Based on 2D hydrodynamic simulations in which the explo- 
sion was triggered by some form of manual energy deposi- 
tion into a stellar progenitor model outside the so-called mass 
cut, they evaluated the composition of the ejecta with a large 
nuclear reaction network. They pointed out that 44 Ti can be 
produced more abundantly in the case o f jet-like exp l osions 
compared to that of spherical explosions. Young et al. (2006) 
examined the composition of the ejecta in three-dimensional 
(3D) SPH simulations to discuss a candidate of the progenitor 
of Cas A. They showed that the abundances of 56 Ni and 44 Ti 
depend on the magnitude and asymmetry of the explosion en- 
ergy as well as on the amount of the fallback. The effects 
of the fallback on the abundances have bee n systematically 
studie d in one-dimensional explosion models (I Young & Frverl 
120071) . More r ecently, 3D effects have been more elab- 
orately studied (Hu ngerford et al.l 2003, 2005), as well as 
the impacts of diff erent explosions by emp loying a num- 
ber of progenitors (Jogg erst et al.ll2009l 1201 Oh or by assum- 
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ing a jet-like explosion dCouch et al.ll2009t lTominagal l2009). 
which is one of the poss i ble candidates of hype rnovae (e.g., 
iMaeda & Nomotol d2003l) : lNagataki et alj d2006l) ). 

In addition to the above-mentioned work, nucleosynthesis 
in a more realistic simulation that models the multidimen- 
sional neut rino-driven SN explos i on has b een also extensively 
studie d dKifonidis et alj I2003L 120061: iGawrvszczak et alj 
2010). Although a small network has ever been included in 
the computatio ns, these 2D simula tions employing a light- 
bulb scheme (Kifonidis et al. 2003) or a more accurate gray 
transport scheme (Sche ck et al.ll2006t iKifonidis et alj|2006l) 
have made it possible to elucidate the nucleosynthesis inside 
from the iron core after the shock-revival up to exp losion in 
a more consistent manner. Kifo nidis et al.1 d2006l) demon- 
strated that the SASI-aided low-mode explosions can most 
naturally explain the masses and distribution of the synthe- 
siz ed elements ob s erved in SN1987A. Their recent 3D results 
bv lHammer etal.1 d2010l) show that the 3D effects that affect 
the velocity of the ejecta as well as the growth of the Rayleigh- 
Taylor instability are really important to correctly determine 
the properties of the ejecta. 

In the present work, we study explosive nucleosynthesis in 
a non-rotating 15 M star with solar metallicity by perform- 
ing 2D hydrodynamic simulations that models a SASI-aided 
delayed explosion via a light-bulb scheme. To extract a de- 
tailed information of the synthesized elements, we follow the 
abundance evolution by employing a large nuclear reaction 
network. It should be emphasized that the mass cut as well as 
the aspherical distribution of the explosion energy are evalu- 
ated from the hydrodynamic simulations, as in our previous 
work on the nucleosy nthesis in magnetohyd rodynamically- 
driven SN explosions dN ishimura et al. 200 g) as well as in 
collapsars dFujimoto et aUl2007l 120081: lOno et alll2009h . 

In §2, this paper opens up with a brief description of a nu- 
merical code for the hydrodynamic calculation, initial condi- 
tions of the progenitor star, and properties of the aspherical 
explosion. In §3, we present a large nuclear reaction network, 
physical properties of SN ejecta, and abundances and masses 
of the ejecta, and heavy-nuclei distribution of the SN ejecta. 
We discuss the uncertainty in the estimate of the abundances 
and masses and compare the evaluated abundances with those 
observed in Cygnus loop in §4. Finally we will summarize 
our results in §5. 

2. HYDRODYNAMIC SIMULATIONS OF AN 
ASPHERICAL NEUTRINO-DRIVEN SUPERNOVA 
EXPLOSION 

2.1 . Hydrodynamic code and initial conditions 

To calculate the structure and evolution of the collapsing 
star, we solve the Newtonian hydrodynamic equations, 

Dp 



D/Dt. The gravitational potential of fluid and the central ob- 
ject with a mass of M; n , $ and $ c , are evaluated with 
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where p,P,\,e, and Y e , are the mass density, the pressure the 
fluid velocity, the internal energy density, and the electron 
fraction, respectively. We denote the Lagrange derivative as 



A$ = 4ttG/9, 



and 



GM n 



(5) 



(6) 



where G is the gravitational constant. We note that M m con- 
tinuously increases due to mass accretion through the inner 
boundary. 

Qe and Qn are the source terms that describe the rate of 
change per unit volume in equations (O and (01, respectively, 
and will be summarized in Appendix A and B. In the present 
study, we take into account absorption of electron and anti- 
electron neutrinos as well as neutrino emission through elec- 
tron and positron captures, electron-positron pair annihilation, 
nucleon-nucleon bremsstrahlung, and plasmon-decays. We 
assume that the fluid is axisymmetric and that neutrinos are 
isotropically emitted from the neutrino spheres with given lu- 
minosities and with the Fermi- Dirac distribution of given tem- 
peratures dOhnishi et al.l2 006). Rates for absorption of neutri- 
nos and neutrino emi ssion through electr on and positron cap- 
tures are taken from Scheck et al. (2006, appendix D). Geo- 
metrical factor f v is set to be 



fu = 



l + Jl-{R v / r f 



(7) 



as in I Scheck et alj (J2006). Here R v is the radius of neutrino 
sphere and is simply estimated with the relation, L v = j^vT* ■ 

4tt R?, for a given set of the luminosity L v and temperature 
T v dOhnishi et alj 120061) . where a is the Stefan-Boltzmann 
constant, We adopt rates for the emission of neutrinos 
(y e ,i> & ,v x ,v^) through pai r annihilation, brem sstrahlung, and 
plasmon-decays as in iRuffert et alj d 19961 Appendix B). 
Moreover, we include the heating term in Qe due to the ab- 
sorption of neutrinos on 4 He and the inelastic scatterings on 
4 He via neutral currents dHaxtonlll98"8l:IOhnishi et alJl2007l) . 

The numerical code for the hydrodynamic calcula- 
tions employed in this paper is based on the ZEUS -2D 
code dStone & Normanl 119921: lOhnishi et al.ll2006l) . We use 
a realistic equation o f state (EOS) ba sed on the relativis- 
tic mean field theory dShen et alj [l998). For lower density 
regime (p < 10 5 g/cm 3 ), where no data is available in the 
EOS table with the Shen EOS, we use another EOS, which 
includes contributions from an ideal gas of nuclei, radiation, 
and el ectrons and positrons with arbitrary degrees of degen- 
eracy (Blinniko v et alJfT996l) . We carefully connect two EOS 
at p = 10 5 g/ cm 3 for physical q uantities to vary cont inuous in 
density at a given temperature (Fujim oto et alj |2006). 

First we perform a spherical symmetric hydrodynamic sim- 
ulation of the core c ollapse of a 15M^ non-ro tating star with 
the solar metalli city dWooslev & Weaverl 1995b using a hydro- 
dynamic code dKotake et al.ll2004l) for aboutlOms after core- 
bounce, when the bounce shock turns into a standing accretion 
shock and the proto neutron star (PNS) grows to ~ 1.2M . 
Then, we map distributions of densities, temperatures, radial 
velocities and electron fractions of the spherical symmetric 
simulation to initial distribution for two-dimensional (2D) hy- 
drodynamic simulations. After the remap, the central region 
inside 50 km in radius i s excised to follow a long-term post- 
bounc e evolution (e.g., IScheck et alJ d2006h : IKifonidis et alJ 
(2006)). We impose velocity perturbations to the unperturbed 
radial velocity in a dipolar manner, and follow the postbounce 
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evolution. The spherical coordinates are used in our simula- 
tions and the computational domain is extended over 50 km < 
r < 50, 000 km and < 9 < tt, or from the Fe core to inner O- 
rich layers, which are covered with 500(r) x 128(6*) meshes. 
The mass is 3.17M Q in the computational domain. We note 
that convective motion occurs at the onset of the 2D simula- 
tion with the above meshes, while the motion does not appear 
in the case of coarser mesh points of 500(r) x 60(6*). Evolu- 
tion of the explosion energy and mass ejection rate are very 
similar to those for high resolution simulation with 300(r) x 
196(61) meshes (50km < r < 3,000km and < 9 < tt) for 
about 350 ms after the core bounce. Therefore resolution of 
the simulations with 500(r) x 128(6*) meshes seems to be ap- 
propriate for the present study. However the resolution may 
be too low to follow later time evolution of the explos ion to- 
wards homologous expansion (Ga wryszczak et aD^OlOh . For 
the high resolution simulation, the minimum grid size in the 
radial- and laterals-directions, Sr and 69, is 1 km and 7r/196, 
respectively, while 5r = 1km and 59 = 7r/128 for our fiducial 
set (i.e., 500(r) x 128(6*) mesh points). 

2.2. Aspherical SN explosion 

We have performed the simulations for models with the 
electron-neutrino luminosities, L Vt =3.7, 3.9, 4.0, 4.2, 4.5, 4.7, 
and 5.0 x 10 52 ergs _1 for 1 -2s after the core bounce, when a 
shock front has reached to a layer with r = 10,000km in al- 
most all directions. We take the input neutrino luminosities 
as above because the revival of the stalled bounce shock oc- 
curs only for models with L Ue >3.9 xl0 52 erg s" 1 , and also 
because for models with L Uf > 5.0 x 10 52 erg s" 1 , the star ex- 
plodes too early for the SASI to grow, as will be discussed 
later. We set L Ce = L Ut and L Vx = 0.5L„ f , where L Ce and L Vx 
are the luminosities of anti-electron neutrino and other-types 
(/i, r, anti-/i, and anti-r), respectively. We consider models 
with neutrino temperatures , T v „, T e „, and T Vr as 4 MeV, 5 MeV, 
and lOMeV, respectively (Ohnishi et al. 2006). The adopted 
neutrino luminosities are comparable to those with a more 
accurat e transport scheme but the temperatures are slightly 
higher dMarek et all 120091; iMarek & .lankul 12009). We will 
present hydrodynamic and nucleosynthetic results for cases 
with lower neutrino temperatures, in §14.11 

We confirm that the explosion are highly aspherical and 
1=1 and / = 2 modes a r e dom i nant as shown i nlKifonidis et al.l 
(120061) : lOhnishi et al.l (120061) : IScheck et al.l 42006J), although 
the shape of the explosion strongly depends on numeri- 
cal d etail, such as mesh resolution and bo undary condi- 
tions dKifonidis et al.l2006tlScheck et al.l2006l) . Entropy con- 
tour of the hydrodynamic simulation is shown in Figure Q]for 
case with L„ e = 4.5 x 10 52 ergs _1 . Most of SN ejecta have 
entropy less than 20£b, where is the Boltzmann constant. 
Entropy attains to 70&b for small amounts of the ejecta. 

We find that for models with L Ve > 3.9 x 10 52 ergs _1 , the 
star explodes aspherically via the neutrino heating aided by 
SASI. Figure [2 a) shows explosion energies as a function of 
L Ue , for all the exploded models. The energies are estimated 
at an epoch of 500 ms after the explosion and slightly increase 
after the epoch. Kinetic and thermal energies of the explosion 
are also shown in Figure|2a). The thermal energies dominate 
over the kinetic ones. 

Higher L Vc makes the onset of explosion earlier and also 
the mass of the PNS (M p ns) smaller. Note that we estimate 
MpNs at f = f exp +500ms, because M p ns only slightly increases 
later than 500 ms after f exp . Here f exp indicates the time scale 
when the explosion sets in, which can be typically estimated 
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FIG. 1. — Entropy contour for case with L Ve = 4.5 X 10" 
after the core bounce. The growth of low-mode SASI is prominent. 
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FIG. 2. — (a) Explosion energies vs. neutrino luminosities. Filled squares, 
circles, and triangles indicate explosion energies, kinetic-part and thermal- 
part of the explosion energies, respectively, (b) M p Ng and f cxp vs. neutrino 
luminosities. Filled squares and circles indicate f exp and M p ns, respectively. 



when the mass ejection rate at 100km grows up O.1M s" 1 in 
our 2D simulations. Figure [2|b) shows the explosion time 
f exp and MpNs as a function of L„ c . Except for the low- 
est luminosity model (L„ f < 3.9 x 10 52 ergs _1 ), M p ns is in 
the range of 1.54- 1.70M Q . These values are much larger 
than the so-ca lled mass-cut in th e sp herical model o f 15Mp 
progenitors by Hashimoto] d 19951) and lRauscher et all d2002l) . 
which is 1.30M Q and 1.32M , respectively. The mass of the 
PNS h owever becomes large r (1 .68M Q ) due to the fallback of 
ejecta dRauscher et al.ll2002l) . 

3. NUCLEOSYNTHESIS IN SUPERNOVA EJECTA 

3.1. Nuclear reaction network and initial composition 

In order to calculate chemical composition of the SN ejecta, 
we need Lagrangian evolution of physical quantities, such as 
density, temperature, a nd, velocity of the material. We adopt a 
tracer particle method (Nagat aki et alJ fl997: Seitenza hl et alJ 
1201 Ol) to calculate the Lagrangian evolution of the physical 
quantities from the Eulerian evolution obtained from our sim- 
ulations. The Lagrangian evolution is followed during the 
2D aspherical simulation as well as the spherical collapsing 
phase. To get information on mass elements, 6,000 tracer par- 
ticles are placed in the regions from 300 extending to 10,000 
km (the O-rich layer). We have confirmed that the estimated 
energies and masses of the ejecta with the 6,000 particles are 
equal to the ones with 3,000 particles within ^1% accuracy 
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and the obtained abundance profiles are also very similar. 

Initial abundances of the particles are set to be those o f 
the star just before the core collapse (Raus cher et al.l l2002), 
in which 1400 nuclei are taken into acc o unt. W e note that 
a presupernova model in Rauscher et al.l (120021) has smaller 
helium, carbon-oxy gen, and oxygen-neon co re masses, com- 
pared with those in Woosl ev" & Weaver! dT995). due to coupled 
effects through the inclusio n of mass loss and the revisions of 
opacity and nuclear inputs (Rausc her et alj|20 02). The mass 
of a particle in a layer is weighted to the mass in the layer. 
We note that the minimum mass of the particles is ~ 10~ 4 M Q . 
We find that more than one fifth particles are ejected due to 
the aspherical explosion. 

Next we calculate abundances and masses of the supernova 
ejecta. Ejecta that is located on the inner region of the star 
('"ej.cc < 10,000km) before the core collapse, has high max- 
imum temperatures enough for elements heavier than C to 
burn explosively. Here r e j cc is the radius of the ejecta at 
the core collapse. We therefore follow abundance evolution 
of the ejecta from the inner region using a nuclear reaction 
net work, which includes 463 nuclide from neutron, proton to 
Kr dFuiimoto et al. 2004). We will discuss effects of neutrino 
interactions on heavy nuclei and uncerta inty in nucle ar reac- 
tion rates on nucleosynthetic results, in {34.2l and l4.31 respec- 
tively. While the abundances of ejecta from the outer region 

(r e j, cc > 10,000km) are set to be those before the core col- and temperatures > 10 1U K. 



The explosion energies also might increase. We emphasize 
that abundances of SN ejecta chiefly depends on the explosion 
energy and the mass of SN ejecta from an inner region, not on 
L Vc , as shown in later. 

3.2. Physical properties of SN ejecta 
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FIG. 3. — (a) Maximum densities vs. maximum temperatures of ejecta 
for model with L Ve = 4.5 X 10 52 ergs -1 . (b) Minimum radial position, 
vs. maximum temperatures of ejecta for model with L„ p = 4.5 X 10 52 ergs -1 . 
Ejecta that falls near the proto neutron star have high densities > 10 9 gem -3 



lapse (Rau scher et al.ll2Q02h . We note that the masses of the 
outer region, or r e j cc > 10,000km, is 10.4M Q . Moreover, 
when temperatures of the ejecta are greater than 9 x 10 9 K, 
we set chemical composition of the ejecta to be that in nu- 
clear statistical equilibrium (NSE), whose abundances are ex- 
pressed with simple analytical expressions, specified by the 
density, temperature and electron fraction. 

Electron fractions of the ejecta are re-evaluated during 
SN explosion coupled with the nuclear reaction network. 
The change in Y e is taken into account through electron 
and positron captures on heavy nuclei, in addition to elec- 
tron and positron captures on neutrons and protons as well 
as absorption of v e and D e on neutrons and protons. The 
captures and absorptions on neutrons and protons are also 
taken into account in the hydrodynamic simulations. The 
rates for the captures and t he abs o rptio ns are a dopted from 
iFuller Fowler. & NewmarJ (Il980i 119821) and IScheck et al l 
(2006), respectively. 

It should be emphasized that post-processing electron frac- 
tions are slightly different (up to 10%) from those estimated 
with hydrodynamic simulations, in which the evolution of 
electron fractions is followed. This is because abundances of 
neutrons and protons in the network calculations are slightly 
different from those estimated with EOS in the hydrodynamic 
simulations. We note that 463 nuclei are taken into account in 
the network calculations, while only neutrons, protons, 4 He, 
and a representative heavier nuclide are evaluated with EOS. 

In the neutrino-heating dominated region, the abundances 
of neutrons and protons in the network calculations are larger 
than those evaluated with EOS. Hence, if we perform hydro- 
dynamic simulations, in which abundances of nucleons are re- 
liably evaluated with the reaction network, the neutrino heat- 
ing rates in the simulations could increase compared to those 
in the current study, since the neutrino heating through the ab- 
sorption of v e and v z is dominant over the other heating reac- 
tions, and the heating rates via the absorption are proportional 
to the abundances of the nucleons. 



Maximum densities, p max , and maximum temperatures, 
r max , are good indicators f or the composition of SN 
ejecta (Thielem ann et al.lll996h . Figure [3 a) shows p max as a 
function of Tmax of the ejecta for L Ve =4.5 x 10 52 ergs _1 . Most 
of the ejecta with relatively low densities (< 10 8 gcm~ 3 ) have 
Pmax and r max similar to those of ejecta in the spherical model 
of core collapse SNe (iThielemann et all l998). For some par- 
ticles that have very high densities simultaneously with high 
temperatures (p max > 10 9 gcm~ 3 and r max > 10 10 K), elec- 
tron captures operate to some extent, so that these particles 
become slightly neutron-rich, T e (10, 000km) < 0.48. Here 
Te(10, 000km) represents the electron fraction for tracer par- 
ticles evaluated when the particles reach r = 10,000km. This 
may be a useful quantity to measure Y e of the ejecta, since Y e 
closely freezes out at r > 10,000km (except for through j3- 
decays at a later epoch). All ejecta with 7; nax > 10 10 K falls 
down to the heating region < 200- 300km to be heated via 
neutrinos (Fig. 0b)). For these neutrino-heated ejecta, p max 
range from 10 8 -2 x 10 10 gcm" 3 . 

For ejecta with higher p max , electron captures on protons 
proceed more efficiently to make its Y e smaller. Time evolu- 
tion of physical quantities of such an ejecta is shown during 
the infall of the ejecta near the cooling region (r < 100km) 
in Figure 0a). As the density and temperature rise to more 
than 10 9 gcm~ 3 and 10 10 K, respectively, the electron fraction 
decreases due to the electron captures. When the ejecta starts 
to be released via neutrino heating at t = 0.42 s, the electron 
fraction increases through the absorption of v e by neutrons. 
Finally Y e (lQ, 000km) becomes 0.461 for the ejecta. 

On the other hand, the electron captures on protons are not 
efficient for a proton-rich ejecta with Y e (lQ, 000km) = 0.559, 
as shown in Figure Hfb). This is because the densities of the 
inner region (r < 200km) are relatively low (< 10 8 gcm~ 3 ) 
due to the mass ejection during an earlier phase (> f exp ). The 
electron fraction therefore remains constant and rises from 0.5 
to 0.559 via the v e absorption in an inner region r < 200km. 
The proton-richness in the ejecta is caused by the small energy 
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FIG. 4. — Time evolution of density, temperature, radial position, 
and electron fraction of an ejecta with (a) y c (10, 000 km) = 0.461 and (b) 
y e (10,000km) = 0.559, for model with h Vt =4.5 x lQ 52 ergs -1 . Solid, dot- 
ted, dashed, and dash-dotted lines represent the density, temperature, radial 
position, and electron fraction of the ejecta, respectively. The electron frac- 
tion is shown with a value multiplied by ten. The density, temperature, and 
the radial position of the ejecta are presented in units of gem -3 , K, and cm, 
respectively. Y e of ejecta can largely change only in an inner region near the 
proto neutron star. 



difference between v e and v e . For T Vc and 7^ adopted in our 
simulations, the r elation, 4(m n -m„) > e Pe — e Vc , holds, which 
leads to Y e > 0.5 dFrohlich et al.l2006al) . where m„ and m p are 
masses of neutron and protons, and e Pc and e„ c are energies 
of anti-electron and electron neutrinos, respectively. We note 
that = 15.8MeV for T Pc = 5MeV and e^ c = 12.6MeV for 
T v . = 4MeV. 




Fig. 5.— 
region r e j, cc 
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Masses as a function of F c (10, 000km) of ejecta from the inner 
< 10,000km for model with L Ve = 4.5 X 10 52 ergs" 1 . 



Figure [5] shows masses as a function of Y e (10, 000km) of 
ejecta from the inner region r e j cc < 10,000km. We find 
that most of the ejecta (98.8%) have electron fractions of 
0.49-0.5. Small fractions of the ejecta, 0.9% and 0.3% in 
mass, are slightly neutron-rich (0.46 <Y e < 0.49) and proton- 
rich (0.5 < Y e < 0.56), respectively. Masses of the slightly 
neutron- and proton-rich ejecta are larger for models with 
larger L v , while the mass fractions of these ejecta are com- 



parable for all the models. 

3.3. Primary and secondary ejecta 
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FIG. 6. — y e (10, 000km) vs. minimum radial position, r m j n , of ejecta for 
model wifhL„ f =4.5 x 10 52 ergs~'. Y e (W, 000km) are largely different from 
0.5 only for ejecta that fall near the proto neutron star. 



Electron fraction of an ejecta with the minimum radial 
position r m i n < 200- 300 km changes due to high neutrino 
flux and/or efficient e ± capture (Fig.Ua)). Figure [6] shows 
F e (10, 000km) as a function of r m ; n of the ejecta for model 
with L Ve = 4.5 x 10 52 ergs _1 . Hereafter, we refer to the ejecta 
with r m ; n < 200km as the primary ejecta, which have high 
maximum densities > 10 8 gcm" 3 and temperatures > 10 10 K 
(Fig.01. The ejecta are heated through the neutrino heating. 
On the other hand, the others are refereed as the secondary 
ejecta, heated chiefly via the shock wave driven by the pri- 
mary ejecta. 

It is true that Y e of the primary ejecta could change if we 
adopt a more accurate neutrino-transfer scheme instead of the 
simplified light-bulb transfer scheme. Abundances of the pri- 
mary ejecta are therefore highly uncertain because of the un- 
certainty on their Y e . On the other hand, for the secondary 
ejecta, Y e changes chiefly through the neutrino absorptions but 
the changes in Y e are found to be less than 1 % for almost all 
the secondary ejecta. In addition, for our typical 2D models 
that produce energetic explosions (L„ f > 4.5 x 10 52 ergs _1 ), 
the masses of the primary ejecta occupy only about 2% 
(8.7 x 1O~ 3 M ) in the ejecta from the inner region (O.41M 
for r e j cc < 10,000km). Therefore, Y e of the secondary ejecta, 
whose mass is much larger than that of the primary ejecta, are 
unlikely to be largely changed even if we use a more accurate 
neutrino-transfer scheme. We conclude that masses of abun- 
dant nuclei, such as lfi O, 28 Si, and 56 Ni, do not largely change 
in the SN ejecta. We will discuss this point in 34.41 

The value of Y e of the primary ejecta depends on the epoch 
of the ejection. The primary ejecta that eject in an early phase 
(before an epoch of 200 -300 ms after the explosion (? e j = 
230 ms)) are mainly neutron-rich, while the primary ejecta 
are proton-rich in the later phase. The proton-rich primary 
ejecta corresponds a kind of neutrino-driven winds, which 
are possibly not neutron-ric h but proton-rich (Fisch er" ;t alJ 
l201ft iHudepohl et alJl20Toh . Figure [7] shows T e (10,000km) 
as a function of f(r max ) of ejecta for model with L v , = 4.5 x 
10 52 ergs _1 , where f(r max ) is defined as the time when the tem- 
perature of ejecta attains to their maximum value, T max . We 
note that the gas starts to be ejected just after t = t(T max ), as 
shown in Figure [4] We note that the dependence of Y e on the 
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FIG. 7. — Y e vs. f(r max ), the time when the temperature of ejecta attains 
to r max for model with L Ve = 4.5 X 10 52 ergs~'. Ejecta at an early phase 
(< 200 — 300 ms after the explosion) tend to be neutron-rich, while proton- 
rich in the later phase. 

epoch of the ejection also appears in a spherical simulation of 
SN explosion of a star with an ONeMg core using an elab- 
orated c ode taking into account an accurate neu trino transfer 
scheme dKitaura et al.ll2006HWanaio et al.ll2009l) . 

3.4. Masses and abundances of ejecta 
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FIG. 8. — Mass of the ejecta from the inne r region and masses of radioac- 
tives vs. L„ e . Note in a spherical model by Rauscher et al. 1 2002) that the 
obtained masses of 56 Ni and 44 Ti are O.IIMq and 1.4 X lO~ 5 A/0, respec- 
tively. 

Masses of nuclei, such as 56 Ni, Ni, Ni and 
44 Ti, have been estimated in some SN remnants; For 
SN1987 A, masses of 56 Ni and 44 Ti are deduc ed to be ~ 
0.07M Q (IShigevama et al.lll988t IWooslevll 19881) and 1-2 x 
10~ 4 M Q (Nagataki 2000, and references therein), respec- 
tively. The estimated mas s of 44 Ti is compar able to that in Cas 
A (1.6!g;| x 1O" 4 M ) dRenaud et al.ll2006l) and greater than 
that in the young est Galactic supernova remnant Gl.9+0.3 
(1-7 x 10 _5 M Q ) dBorkowski et al.ll2010l) . which may origi- 
nate from a Type la event. Figure[8]shows masses of SN ejecta 
ejected from the inner region (r e j, cc < 10,000km), M e j i n , and 
masses of 56 Ni and 44 Ti. Masses of the ejecta from the inner 
region of r e j cc < 10, 000km (solid line with filled squares) and 
masses of 56 Ni (dashed line with filled circles) and 44 Ti (dot- 
ted line with filled triangles) are shown with a value times a 
factor of 1, 10 and 10 4 , respectively. We find that masses of 
the ejecta and 56 Ni roughly correlate with the neutrino lumi- 



nosities. The masses of Ni and Ti a re less than and com - 
parable to those in the spherical model (Raus cher et al.ll2002l) . 
or 0.1 IM Q and 1.4 x 10" 5 M Q , respectively. Note that the ex- 
plosion energy is 1 .2 x 10 51 erg in the model. We find that ^Ti 
relative to 56 Ni are much smaller than those in the solar sys- 
tem, in SN1987A, and in Cas A, but comparable to that in the 
supernova remnant Gl.9+0.3. 57 Ni relative to 56 Ni is compa- 
rable to those in the solar system and in SN1987A, while 58 Ni 
relative to 56 Ni is overproduced compared with those in the 
solar system and in SN1987A, because the y are abundantly 
produced in the sligh tly neutron-rich ejecta ( Hashimot o! 1 9951; 
lNagatakietaDl997l) . 

It should be emphasized that 44 Ti is underproduced in our 
simulations of the SN explosi on, contrary to the o verproduc- 
tion in the previous 2D results (Nag ataki et al.l 19971) . in which 
the explosion energy is aspherically and artificially added and 
the remnant mass is set to a value in order that the ejected 
m ass of 56 Ni is reprodu ced the mass observed in SN1987A 
in Nagatak i et al.1 (119971) . The underproduction is possibly 
caused by a larger remnant mass in our models. This is be- 
cause the explosion energy and the ratio of the explosion en- 
ergy on the polar axis to that on the equatorial plane are com- 
parable to those evaluated in our simulation, and masses of 

Ti as well as 56 Ni have shown to strongly de pend on the 
value of the remnant mass in 2D calculations (lYoung et al.l 
2006). 
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FIG. 9. — Overproduction factors, X//Xqj vs. mass number for (a) L„ t , = 
4.5 X 10 52 ergs -1 and for (b) L„ e = 3.9 X 10 52 ergs -1 Thick horizontal-dashed 
lines represents a factor equals to that of 16 0, while two normal and two thin 
lines denote a factor equals to that of 16 times 2, 1/2, 3, and 1/3, respectively. 

In order to compare estimated abund ances with the solar 
system ones (lAnders & Grevessdfl989l) . we have integrated 
masses of nuclei over all the ejecta to evaluate the abun- 
dances of the SN ejecta. Figure |9j a) shows overproduction 
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factors after decays as a function of the mass number, A for 
L Ue = 4.5 x 10 52 ergs _1 , in which the explosion energy is 1.3 
xlO^ergs" 1 (Fig.Eta)) and M ej , in is O.41M (Fig.©. We 
find that the abundance pattern of the SN ejecta is similar to 
the solar system one. We note that 17 0, which is underpro- 
duced in the ejecta, can be abundantly synthesized in Type 
la SNe, and that 15 N and 19 F are com parably produced if 
neutr ino effects are taken into account (IWoosley & Weaver! 
1995), as shown later in Figure [14] We poin t out that 64 Zn, 
which is underproduced in the spherical case (Rausch er et alj 
2002), is abundantly produced in slightly neutron-rich ejecta 
(0.46 < Y e < 0.49). 

In the ejecta, however, the neutron-rich Ni, 62 Ni, is over- 
produced. The overproduction of the neutron-r ich Ni iso- 
topes has also appeared in t he spherical models (Hashimoto 
Il995t iRauscher et all 12 002) as well as in the asymmetric 
model (Nagata ki et al.lll997h . The overproduction has been 
suppressed if the electron fraction of slightly neutron-rich 
material (Y e = 0.495 just before the core collapse) near the 
mass cut is artificial ly modified to 0.499 (Hashimoto] 1 1 995b 
Nagat aki et alJI 1997b . Hence, it has been remarked that the 
overabundance of the neutron-rich Ni isotopes may inherit the 
uncertainty in the progenitor model, in particular in the neu- 
tronization due to the combined effects of convective mixing, 
electro n captures, and p ositron decays during their Si burring 
stage ( Hashimotol l 1 9951) . In our aspherical models, the over- 
abundances of Ni in the primary neutron-rich ejecta are not 
large, other than 60 Ni, as we will discuss later Q4.4I ). The 
overproduction of Ni mainly takes place in slightly neutron- 
rich secondary ejecta (0.49 < Y e < 0.5). We find that the 
electron fractions of the secondary ejecta are 0.4985 just be- 
fore the core collapse and decrease by up to 1% through 
the neutrino absorptions. Therefore, the overabundance of 
62 Ni possibly inherits the uncertainty not only in the pro- 
genitor model but also in the change for Y e during the SN 
explosion. Moreover, the overproduction of Ni isotopes are 
also shown to be reduced in the ejecta from the spherical 
SN explosion induced via artificially enhanced neutrino heat- 
ing dFrohhch et alj|2006al) . 

On the other hand, the integrated abundances are largely 
different from those in the solar system for L Ve = 3.9 x 
10 52 ergs _1 , in which E exp = 0.45 x 10 51 ergs _1 and M e ^ in = 
0.20M Q . Nuclei with A > 28 are deficient in the model 
(Fig.HJb)), compared with O, which is mainly ejected from 
outer layers. The overproduction factors of nuclei lighter than 
27 Al are very similar for both models, because these nuclei 
are mainly synthesized during the hydrostatic burning. While 
nuclei heavier than Si are mainly produced during the SN ex- 
plosion through the explosive burning in ejecta from the inner 
region (r e j, cc < 10,000km). 

The abundance patterns for models with M e j ; n ~ (0.4- 
0.5)M o , or those with L Ve = (4.5, 4.7, and 5.0) x l0 52 ergs _1 , 
are similar to that in the solar system. We should notice 
that M e j.in, which corresponds to masses of a neutron star, 
estimated as (1.54- 1.62)M for these models (Fig.^b)), 
are comparable to the baryonic mass (around 1.5M Q ) of a 
neutron star for observed neutron-star binaries, in which the 
progenit or of the neutron s tar seems to be a star with a 
Fe core dSchwab et alj|2010h . Moreover, E exp /M ei = (0.92- 
1.2) x 10 50 erg/M Q evaluated for the models are comparable 
to the v alue E exp /M e ; = 0.76 x 10 50 e rg/M n estimated in SN 
1987A dShigeyama & Nomotol 19901) . The time from the core 
collapse to the explosion, f exp , is ~ 0.2 s (Fig.|2fb)), which is 



enough to grow a low-mode SASI. The growth is appropriate 
for explaining the distributions and velocities of nuc lei ob- 
served in SN 1987A, as shown in lKifonidis et alj d2006l) . 

3.5. Aspherical infall 
The value of M e j ; n depends on E exp as well as aspherical 
matter infall to the neutron star, while the masses of 56 Ni cor- 
relates with £exp (Figs.|2|a) and©. Figure [10] shows the posi- 
tion of SN ejecta before the core collapse for model with L v , = 
4.5 x 10 52 ergs _I . Filled squares, filled circles, filled triangles, 
and open triangles indicate ejecta, which correspond to the 
complete Si burning (T mdXt9 > 5 where r max ,9 = 7m ax /10 9 K), 
the incomplete Si burning (4 < r max 9 < 5), the O burning 
(3.3 < r max , 9 < 4), and the C/Ne burning (1 < r raax . 9 < 3.3), 
respectively. We find that whole of the iron core (~ 2,000 km) 
collapses to the proto neutron star and that the infall of mate- 
rial to the proto neutron star is aspherical; Larger amounts of 
the gas infall from the upper hemisphere compared with from 
the lower half. Larger amounts of material aspherically infall 
to the neutron star for model with L„ f = 4.5 x 10 52 erg s" 1 com- 
pared with model with L„ f = 5.0 x 10 52 ergs _1 . Though £ exp is 
comparable, M e ji n for L Vf = 4.5 x 10 52 ergs _1 is smaller than 
that for L Ue =5.0 x lO^ergs" 1 . 

1-0 awr 1 1 
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FIG. 10. — Initial position of SN ejecta that are released from the inner 
region (r cc < 10,000km) for model wifhL„ t , =4.5 X 10 52 ergs~'. 



3.6. Aspherical distribution of energy and nuclei 

Distributions of the energy and abundances of SN ejecta are 
highly aspherical, as one can expect from the entropy distri- 
bution (Fig.[T|l. The distributions are shown for L„ e = 4.5 x 
10 52 ergs _1 at t = 1.5s in Figure[TT| The energy and abun- 
dances are presented in units of 10 47 erg and normalized by the 
solar abundances, respectively, for logarithm scale. It should 
be emphasized th at the distribu t ions change during later 
expan sion phase dKifonidis et al.1 l2006t iGawrvszczak et all 
2010). We note that the secondary ejecta that is located in 
outer layers (r e j cc > 10,000km) are not shown in the fig- 
ures. The high energy ejecta concentrates on the shock front, 
in particular in polar regions (9 < ir/4,6 > 3ir/4). Both Si 
and O burning proceed to produce 56 Ni and 28 Si abundantly 
in the region. We emphasize that a deformed shell struc- 
ture forms with Fe, Si and O layers from inside to outside 
in the secondary ejecta heated via the shock wave. Compo- 
sition of the secondary ejecta depends chiefly on T m - dx , which 
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FIG. 11. — Distribution of the energy and abundances of the ejecta for L Ve = 4.5 X 10 52 ergs 1 at t = 1.5 s in a region of < x < 20,000km and -20,000km < 
y < 20,000km. Abundances are shown in units of those in the solar system for Ifi O, 28 Si, and Fe, which is the sum of 56 Fe, 56 Ni, and 54 Fe. 



is determined via the shock heating. The shell structure is 
therefore formed in the secondary ejecta, because T mdx gradu- 
ally decreases from inner to outer layers, where the structure 
is highly aspherical and deformed. We note that Ca is co- 
produced with 28 Si, as in the spherical models. 

The primary ejecta however do not form such a structure. 
This is because physical properties of the ejecta, which are 
mainly heated via neutrinos to be > 10 10 K (Fig. 0J, are in- 
dependent from their positions, and thus their compositions 
are also independent. A fraction of the primary ejecta that 
are mainly composed of Fe and Ni is found to be mixed into 
the deformed layers, while a small amount of the secondary 
ejecta falls towards the neutron star through near the equato- 
rial plane. Consequently, the ejecta with abundant O and/or 
Si appear in the central region (r < 2,000km). Due to such 
mixings and infall, some of the Si- and O-rich ejecta can pen- 
etrate deeper into the Fe-rich ejecta, which is in sharp contrast 
to spherical models. 

We have averaged energy and masses of nuclei of the in- 
ner ejecta over a radial direction. We note that the energy 
and masses of the outer ejecta (r e j cc > 10,000km) are not in- 
cluded in an averaging procedure. If we include ejecta from 
the outer layers, the asymmetry of lfi O becomes small, be- 
cause of the existence of spherically distributed 16 in the 
layers. Figure \l2\ a) shows the averaged energy and masses 
of lfi O, 28 Si, and 56 Ni. Asphericity of the energy is promi- 
nent; Ratio of the energy at 9 = tt/2 to that at 9 = is ^4 
and ratio of the maximum to the minimum energy is ~ 8. We 
note that the ratios are compar able to those in two-dimension 
models dNagataki et al .11 19971 models A2 and A3). Only a 
small amount of 56 Ni exists near the equatorial plane (9 ~ 
(0.8 - 1 .6)rad.), where the energy is lower compared with that 
near the polar axis. The asymmetry of 28 Si mass is similar to 
that of the energy and smaller than that of 56 Ni. Mass distribu- 
tion of 44 Ti is very similar to that of 56 Ni, as in spherical mod- 



els, because both nuclei are produced in the secondary ejecta 
through a-rich freezeout. The peaks of 28 Si mass around 
the polar directions are misaligned with those of 56 Ni around 
9 ~ 7r/4 and 37r/4, although the misalignment could disap- 
pear due to the lateral motion of the ejec ta, during later ex- 
plosio n phase (< 300 s) as pointed out in iGawrvszczak et al.l 
(2010). We note that the lateral motion does not follow during 
the later phase in the present study. 

In addition to the mass fractions of nuclei, number ratios of 
nuclei relative to O ( ie O) are important quantities, because the 
ratios are indicative to the production and destruction mech- 
anism of nuclei. Angular distributions of the number ratio 
are shown for the inner ejecta (r e j cc < 10,000km) averaged 
over a radial direction for L Ve = 4.5 x 10 52 ergs _I in Figure 
[T2rb). Ratio of Ne ( 20 Ne, solid line) is 0. 1-0.2 for a direction 
where the O/Ne burning does not effectively operate, while it 
becomes much lower in a region where the O/Ne burning ef- 
fectively operates to synthesize abundant Si and/or Fe relative 
to O. Distribution of the number ratio N(Mg)/N(0) ( 24 Mg, 
dash-dotted line) is similar to the spherical case, because Mg 
is produced during both the stellar evolution and the explosion 
through the same process, or the O/Ne burning, with which 
2V(Mg)/iV(0) becomes ~ 0.05 for a 15M Q progenitor. Here 
N(X) is the number of nuclei, X. Ratio of Si ( 28 Si, dashed 
line) anti-correlates with that of Fe (the sum of 56 Ni, 56 Fe, 
and 54 Fe, dotted line). This is because Fe is produced via the 
burning of Si. 

4. DISCUSSIONS 

4. 1 . Dependences on neutrino temperatures 

We have assumed the neutrino spheres with given lumi- 
nosities and with the Fermi-Dirac distribution of neutrino 
temperatures, and considered models with neutrino temper- 
atures, T v , T(j, and T Vx as 4MeV and 5MeV, and lOMeV, 
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FIG. 12. — (a) Angular distribution of the averaged energy and masses 
of nuclei of the inner ejecta (r e j cc < 10,000 km) over a radial direction for 

L Ve = 4.5 X 10 52 ergs~' at t = 1.5s. Solid, dashed, dotted, and dash-dotted 
lines indicate the averaged mass of 16 0, 28 Si, and 56 Ni and the averaged 
energy, respectively. The masses and energy are shown in units of O.O1A/0 
and 10 erg, respectively, (b) Angular distribution of number ratios of nuclei 
to I6 Oofthe inner ejecta (r c j cc < 10, 000 km) averaged over a radial direction 
for L„ c = 4.5 X 10 52 erg s _1 at t = 1 .5 s. Solid, dashed, dotted, and dash-dotted 
lines indicate the ratios of 20 Ne, 28 Si, iron, and 24 Mg, respectively. The iron 
consists from 56 Ni, 56 Fe, and 54 Fe. 



model, in spite of only 10% changes in neutrino temparatures. 
Table 1 summarizes hydrodynamic and nucleosynthetic prop- 
erties for the nine models that produce explosions, in which 
T v s are changed in the three ways as mentioned above. The 
explosion energies are (0.35, 0.35, and 0.45) xl0 51 erg for 
models with L Ve = (4.0, 4.5, and 5.0) xl0 52 ergs _1 , respec- 
tively. Mass of the ejecta from the inner region (r e j cc < 
10,000km), M e j, in , and 56 Ni mass are O.11-O.21M and 
0.01 1 -0.017M Q , respectively. Figure Q~3] shows overproduc- 
tion factors for the low T v model with L Ve = 5.0 x 10 52 ergs _1 . 
We emphasize that the overproduction factors are very similar 
to those for the high T v model with L Ve = 3.9 x 10 52 ergs _1 , in 
which the explosion energy and the ejecta mass from the inner 
region are comparable to the above model, or 0.45 x 10 51 erg 
and 0.020M Q , respectively. We conclude that the abundances 
of the ejecta do not directly depend on T v but mainly depend 
on the explosion energy and the ejecta mass. 
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FIG. 13. — Same as Figure|9]but for model with low neutrino temperatures 
andL„ e = 5.0 X 10 52 ergs _1 . The factors are very similar to those for a model 



with high neutrino temperatures and L„ f 
Figure[§b). 
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TABLE 1 

Model parameters and properties of the explosion for nine 
exploded models 





To. 






*exp 


MpNS 


Mej iin 


M( 56 Ni) 


M( 44 Ti) 


4.0 


5.0 


3.9 


0.45 


1.20 


1.84 


0.20 


1.85e-2 


5.29e-6 


4.0 


5.0 


4.0 


0.80 


0.31 


1.70 


0.34 


2.72e-2 


1.31e-5 


4.0 


5.0 


4.2 


1.02 


0.31 


1.65 


0.39 


4.51e-2 


8.04e-6 


4.0 


5.0 


4.5 


1.30 


0.23 


1.62 


0.41 


5.49e-2 


1.52e-5 


4.0 


5.0 


4.7 


1.03 


0.21 


1.54 


0.49 


5.10e-2 


1.14e-5 


4.0 


5.0 


5.0 


1.30 


0.18 


1.56 


0.48 


5.44e-2 


1.05e-5 


3.6 


4.5 


4.0 


0.35 


1.41 


1.93 


0.11 


1.10e-2 


6.95e-6 


3.6 


4.5 


4.5 


0.35 


1.19 


1.87 


0.17 


1.40e-2 


5.53e-6 


3.6 


4.5 


5.0 


0.45 


0.62 


1.83 


0.21 


1.73e-2 


6.38e-6 


Note. — 


Each column shows 7^ , 


j 1 - t l u 


^expi ^exp 


and masses 


(M pNS ,Mq,i„ 



M( 56 Ni), and M( 44 Ti)), in units of MeV, MeV, 10 52 ergs"', 10 51 erg, s, 
respectively. 



respectively, The temperatures are however slightly lower 
and change in time, as s hown in simulations using a more 
elabo rate numerical code (M arek et al.ll2009i: iMarek & Jankal 
2009). In order to evaluate dependences of hydrodynamic 
and nucleosynthetic results on neutrino temperatures, we have 
performed simulations for three models with lower neutrino 
temperatures; T Ue = 3.6 MeV, T$ e =4.5 MeV, and 7^ = 9MeV, 
andL^ = (4.0, 4.5, and 5.0) Xl0 52 ergs _1 . We find that, for 
a given L Ue , the explosions are much weaker for a lower T v 



4.2. Effects of the v interactions on heavy nuclei 

Effects of neutrino interactions on S N ejecta have been in- 
vestigated for spherical S N explosion ( Wo oslev et al.l [19901: 
IWooslev & WeaveJ ll995). We have calculated abundances 
of the ejecta taking i nto ac count v interactions on heavy 
nuclei fcoriely et all 1200 ll) for L Ve = 4.5 x 10 52 ergs _1 . 
Cross section s for t he neutrino reactions are taken from 
IWooslev et al.l dl990l) . in which rates are included for neu- 
tral and charged current reactions on He and from C to Kr. 
Figure [14] shows ratios of the abundances of the ejecta with 
the v interactions to those without the interactions. We 
find that the increase of the abundances due to the interac- 
tions of v on heavy nuclei are up to 50%, except for 19 F, 
whose ratio is 2.1. 19 F is enhanced via the neutral and 
charged current interactions of 20 Ne, or 20 Ne(z/z; ,p) 19 F and 
20 Ne(f'e,e + n) 19 F, respectively. 15 N and 50 V are synthesized 
through I6 0(w',p) 15 N and 50 Cr(£> e ,e + ) 50 V, respectively. We 
note that the enhancement of light elements, such as Li, Be, 
and B, via the v interactions is not fully taken into account in 
our calculation, because the composition of ejecta from outer 
layers (r e j cc > 10,000km) is fixed to be the pre-SN composi- 
tion, although the production of these element s is efficient in 
the outer layers through n eutrino interactions (IWooslev et al.l 
Il990t Yoshida et al. 2008). Moreover, neutrino absorptions 



10 



Fujimoto et al. 



2.2 |— 

2.0 

1.8 - 

g 1.6 - 

>? 1.4 - 

1.2 

1.0 ■ 

0.8 — 
10 



ft? 



30 40 50 

Mass number 



60 



FIG. 14. — Ratios of the abundances of SN ejecta taking into account v 
interactions on heavy nuclei to those without the u interactions, as a function 
of A for L Ve = 4.5 X 10 52 erg s -1 . 

on D may be important for n ucleosynthesis as well as the dy- 
namics of the SN explosion (Nakamura et al. 20Q9|). 

4.3. Dependences of abundances on nuclear reaction rates 

We have adopted reaction rates mainly taken from the REA- 
CLIB database in our nuclear reaction network presented in 
33.11 The database is recently updated and continuously 
maintained by the Joint I nstitute for Nuclea r Astrophysics 
(JINA) REACLIB project dCvburt et alj|2010t) . We have cal- 
culated abundances of the ejecta for L Vt = 4.5 x 10 52 ergs _1 
adopting reaction rates taken from the JINA REACLIB VI. 
database. Figure Q3] shows ratios of the abundances of the 
ejecta adopting reaction rates in the JINA REACLIB V1.0 
database to those in REACLIB database. We find that the 
differences are small, up to 30%, between abundances with 
the JINA REACLIB V1.0 and REACLIB databases. If we 
use newly evaluated q-captur e rates on 40 Ca and 4 4 Ti, yields 
of 44 Ti are likely to be lower (iHoffman et alJlloToh . 
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FIG. 15. — Ratios of the abundances of SN ejecta using reaction rates in the 
JINA REACLIB VI. database to those in REACLIB database, as a function 
of A for L Ue = 4.5 X 10 52 erg s" 1 . 



4.4. Uncertainty of abundances of the ejecta 

As shown in 33.31 the ejecta consists of the primary ejecta, 
whose Y e and thus abundances are highly uncertain and the 
secondary ejecta, which have a relatively definite composi- 
tion and are much heavier than the primary ejecta. In or- 
der to clarify uncertainty on the estimate of abundances of 



the ejecta, we have evaluated abundances of the ejecta be- 
ing integrated over ejecta without the primary ejecta and have 
compared above the abundances with those summed over all 
the ejecta (Fig.[9ja)). Figure [16] shows ratios of the abun- 
dances of all the ejecta to those without the primary ejecta 
for L Ue = 4.5 x 10 52 ergs _1 . We find that abundances of most 
of nuclei with A < 70 do not largely change within a factor 
of 2. The overproduction factors are therefore very similar 
to those in Figure |3 a )- The ratios are greater than 1.5 for 
43 Ca, 45 Sc, 47 Ti, 50 Ti, 54 Cr, 60 Ni, 64 Zn, 66 Zn, and ™Ge, in 
particular, 9.0, 2.8, and 3.1 for 64 Zn, 66 Zn, and 70 Ge, respec- 
tively. These nuclei are therefore chiefly synthesized in the 
primary ejecta, not in the secondary ejecta. Moreover, 50 Ti, 
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FIG. 16. — Ratios of the abundances of SN ejecta to those without the 
primary ejecta, as a function of A for L Ue = 4.5 X 10 52 erg s _1 . 

54 Cr, 60 Ni, 64 Zn, 66 Zn, and 70 Ge are synthesized in neutron- 
rich primary ejecta, while 43 Ca, 45 Sc, 47 Ti, and 60 Ni are abun- 
dantly pro duced in proton-rich pr imary ejecta through vp- 
processes dFrohUchetal J 1200635 IPruet et al J 120061: [Wanaio 
l2006h . 

In short, abundances of 64 Zn, 66 Zn, and 70 Ge are highly 
uncertain but those of the other nuclei are relatively definite. 
If a fraction of the primary ejecta could become much larger, 
abundances might be highly uncertain, in particular for the 
nuclei that are produced in the primary ejecta. The fraction is 
however unlikely to be much larger, because the ejection of 
the secondary ejecta is driven by the shock wave caused by 
the primary ejecta. It should be noted that the mass fractions 
of the primary ejecta to the secondary ejecta are comparable 
for all the models, although the explosion energies and ejected 
masses from the inner region diverse among the models. 

Moreover, matter near the proto neutron star is blown off 
via neutrino-driven winds during later evolution of the star 
(> 2s). The ejecta could be not neutron-rich b ut proton-rich 
(Y e ~ 0.5 - 0.6) and have hig h entropy (> 50k B ) dFischer et all 
l2010tlHudepohl et al.ll2010h . The z^p- process could operate in 
the ejecta to synthesize light p-n uclei dFrohlich et alj2006albt 
IPruet et al J 120061: IWanaiol 120061) . However, the process may 
have not large contributions to the ejected masses of abundant 
nuclei, since the mass ejection rate is sm all through the winds 
at the later epoch (Hiidepohl et al.ll2010t) . 



4.5. 



Comparison with observations on abundances of SN 
remnants 



Recently, iKimural d2009t) : lUchida et al] d2009l) have ana- 
lyzed the metal distribution of the Cygnus loop by using the 
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data obtained by the Suzaku and XMM-Newton observations. 
It is pointed out that the progenitor of the Cygnus loop is 
a core-collapse supernova explosion whose progenitor mass 
ranges in ~ 12- 15M Q . The ejecta distributions are asymmet- 
ric to the geometric center: the ejecta of O, Ne are distributed 
more in the north-west rim, while the ejecta of Si and Fe are 
distributed more in the south-west of the Cygnus loop. Since 
the material in this middle-aged supernova remnant has not 
been completely mixed yet, the observed asymmetry is con- 
sidered to still remain a trace of inhomogeneity produced at 
the moment of explosion. These evidences may allow us to 
speculate that the asymmetry of the heavy element observed 
in the Cygnus loop may come from globally asymmetric ex- 
plosions explored in this work. We try to seek the relevance 
in the following. 

Figure [17] shows number ratios of Ne, Mg, Si, and Fe rel- 
ative to O observed in Cygnus loop and those of ejecta for 
L Ue = 4.5 x 10 52 ergs _1 . We find that the ratios observed in 
Cygnus loop (solid line) are comparable to those averaged 
over ejecta from the inner lower region (r e j cc < 10,000km 
and 9 > 2ir/3) (dotted line). The ratios averaged over all 
the ejecta (dash-dotted line) are same as those averaged over 
ejecta from the inner equatorial region (r e j cc < 10, 000 km and 
it/3 < 9 < 2"7r/3) (dashed line). We note that the ratios av- 
eraged over the ejecta from the inner upper region (r e j cc < 
10,000km and 9 < tt/3) are comparable to but slightly lower 
than those averaged over the ejecta from the inner lower re- 
gion (dotted line). 
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FIG. 17. — Number ratios of Ne, Mg, Si, and Fe relative to O observed 
in Cygnus loop and those of ejecta for L Ve = 4.5 X 10 52 ergs~' at t = 1.5 s. 
Solid line indicates ratios observed in Cygnus loop, while dotted, dashed, and 
dash-dotted line indicate ratios averaged over ejecta from an inner equatorial 
region (r e j cc < 10.000km and n/3 < 9 < 2?r/3), ejecta from an inner polar 
region (r c j cc < 10,000km and 9 > 27r/3), and all the ejecta, respectively. 

7Y(Ne) /N(O) and N(Mg) /N(O) are independent from an av- 
eraging region and are comparable to all our aspherical mod- 
els. These ratios are therefore concluded to be determined 
during the hydrostatic evolution of the progenitor and thus are 
a good indicator for the progenitor mass. On the other hand, 
ratios of Si and Fe relative to O are much higher in the inner 
region, in particular in the inner lower region, compared with 
the outer region (> 10,000 km). N(Si)/N(0) and /y(Fe)//Y(0) 
averaged over all the ejecta are much smaller than those ob- 
served in Cygnus loop. Therefore incomplete radial mixing of 
Fe and Si with O is required to be explained the high Fe and Si 
ratios in Cygnus loop. In fact, the averaged abundances of the 
Cygnus loop have a correlation with the Si- and Fe- rich re- 



gions dKim ura 2009; Uchida et alj |2009l). Our results suggest 
that the abundance ratios for Ne and Mg as well as Si and Fe 
observed in the Cygnus loop could be well reproduced with 
the SN ejecta from the inner region of a 15M progenitor. 

4.6. Three dimensional effects 

In order to draw a robust conclusion to the find- 
ings obtained in the current 2D simulations, it is indis- 
pensable to move on to three dimensio n al (3D) simu- 
lations (e.g .Blondin & Mezzacappa (2007); Iwaka mi et al. 
(120081 120091) : iNordhaus et al.1 (2010); Wongwathanarat et al. 
(120101) 1. In 2D, the growth of SASI and the large-scale con- 
vection tend to develop along the coordinate symmetry axis 
preferentially, thus suppressing the anisotropics in explosions 
as well as in the resulting explosive nucleosynthesis. An en- 
couraging news to us is that 3D simulations cited above are, 
at least, in favour of the SASI-aided low-modes explosions, in 
which the major axis of the expanding shock has a certain di- 
rection at the moment of explosion (albeit free from the coor- 
dinate symmetry axis). By incorporating the presen t scheme 
to the 3D simulations of llwakami et all d2008ll2009h . we plan 
to clarify these issues as a sequel of this study. 

5. SUMMARY 

We have investigated the explosive nucleosynthesis in the 
delayed neutrino-driven, aspherical SN explosion aided by 
SASI, based on 2D, axisymmetric hydrodynamic simulations 
of the explosion of a non-rotating 15M star. We employed 
a hydrodynamic code with a simplified light-bulb neutrino 
transport scheme. We have approximately taken into account 
neutrino heating and cooling as well as the evolution of elec- 
tron fraction due to weak interactions, both in the hydrody- 
namic simulations and nucleosynthetic calculations. Neutri- 
nos are assumed to be isotropically emitted from the neutrino 
spheres with given luminosities and with the Fermi-Dirac dis- 
tribution of given temperatures. We have performed simula- 
tions with the temperatures and luminosities of v e , v e and v x 
constant in time. We have followed abundance evolution of 
SN ejecta using the nuclear reaction network coupled with an 
evolution equation of the electron fraction of the ejecta. 

We summarize our results as follows, 

1. The stalled shock revives due to the neutrino heating 
aided by SASI for cases with L„ e > 3.9 x 10 52 ergs _1 
and the aspherical shock passes through the outer lay- 
ers of the star (> 10,000km). Evaluated explosion en- 
ergies roughly correlate with neutrino luminosities. For 
models with larger luminosities, the explosion occurs 
earlier and mass of a neutron star becomes lighter. 

2. Whole of the iron core of the progenitor collapses to the 
proto neutron star. The infall of material to the star is 
aspherical. Larger amounts of the gas infall from direc- 
tions with lower explosion energies. 

3. Abundances of the neutrino-heated ejecta are highly un- 
certain, in particular for neutron-rich ones, due to the 
uncertainty in the estimate on Y e . On the other hand, 
the shock-heated ejecta has definite abundances, which 
depend on the maximum temperature mainly. The un- 
certainty in the estimate of the masses and abundances 
of abundant nuclei in the SN ejecta is small because of 
the small fraction of the neutrino-heated ejecta. 
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4. Abundance pattern of the supernova ejecta is similar 
to that of the solar system, for cases with the mass of 
the ejecta from the inner region (< 10, 000km), M e j i n = 
(0.4-0. 5)M©, which corresponds to models with a high 
explosion energy of ~ 10 51 erg. Masses of a neutron 
star remnant, estimated to be (1 .54- 1 .62)M Q for these 
models, are comparable to the baryonic mass of the 
neutron star observed in neutron-star binaries. E exp /M e j 
evaluated for the models ((0.92-1.2) x lO 5O erg/M ) 
are comparable to the estimate in SN 1987 A. 

5. Underproduction of 44 Ti and overproduction of 62 Ni, 
which appear in spherical models, are also shown in our 
2D calculations. The overproduction of 62 Ni possibly 
inherits the uncertainty not only in the progenitor model 
but also of the change in Y e during the SN explosion. 
On the other hand, 64 Zn, which is underproduced in a 
spherical model, is found to be abundantly produced in 
our 2D model, although the abundance and mass are 
uncertain. 



6. Distributions of nuclei and energy are highly aspherical 
in the SN ejecta, although the progenitor is non-rotating 
and has spherical symmetric configuration. The shock- 
heated ejecta forms aspherical and deformed shell-like 
structure composed of Fe, Si, and O from inside to out- 
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nite structure and fractions of the ejecta are mixed into 
the shell-like structure of the shock-heated ejecta. 



7. The asymmetry of Si mass distribution is similar to 
that of the explosion energy and smaller than that of 
56 Ni. 40 Ca and 44 Ti are accompanied with 28 Si and 56 Ni, 
respectively. 



8. The ratios for Ne, Mg, Si and Fe observed in Cygnus 
loop are well reproduced with the SN ejecta from the 
inner region of the 15M progenitor. The number ratios 
for Ne and Mg relative to O seem to be good indicators 
for the mass of a progenitor of a core-collapse SN. 
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APPENDIX 

In this appendix, we summarize our treatment on the rate of change in Y e and energy per unit volume, £>n and Qe, respectively. 
We take into account absorption of electron and anti-electron neutrinos as well as neutrino emission through electron and positron 
captures, electron-positron pair annihilation, nucleon-nucleon bremsstrahlung, and plasmon-decays. Moreover, we incl ude the 
heating term in Qe due to the absorption of neutrinos on 4 He and the inelastic sca tterings on 4 He via neutral currents (iHaxtonl 
' 1988t lOhnishi et al.ll2007l) . We chiefly follow the treatments in Appendix D of [Scheck et alJ (120061) and in Appendix B of 
Ruff ertet all (1 19961) and references therein). 

THE RATE OF CHANGE IN SPECIFIC ENERGY, Q E 

neutrino absorption processes 

The heating rate per unit volume through the absorption of v t on neutrons is described as 

Ue = ac ^K (o e((e - }) ' (A1) 

where A = (m n -m p )c 2 , n„ is the number density of neutrons, and a = \G\rn\h 2 /ire 2 = |(3a^,+ l)ao/(m e c 2 ) 2 , with the Fermi 
coupling constant Gf, the reduced planck constant Ti, the electron mass m e , a w = 1 .254, and ctq = 1 .76 x 10~ 44 cm 2 . Here, the n-th 
energy moment of z/ e , (e" ), is given by 

«) = (W"^^, (A2) 
and a factor for the Pauli blocking, 0( (£„,.)), is approximated as 

e((e, j ))=l-/FD(-^^,r ?e -Y (A3) 

where r\ Vl is the chemical potential of vi in units of k^T and set to be 0, ry e - is the chemical potential of electrons in units of k^T, 
ks is the Boltzmann constant, and J-„{rj) is defined as, 



T n {ri)= / dxx"/ FD (x,?7) (A4) 
Jo 

using the Fermi-Dirac distribution function, 

fm(x, v) = t- ; r ■ (A5) 

1 + exp(x-?7) 

The heating rate via the absorption of D e on protons is given by 

L PeWp (4 c r+3A(4r + 3A 2 ( £Pc r+A 3 ( e g e r 
e = ac ^cj- c K) ' tA6) 

where the n-th energy moment of D e , (ej.)*, is given by 



Capture processes 

The energy emission rate of v e per unit volume via electron capture on protons is given by 



Ql c = — n p n e - [(eiy + 2A{el-y + A 2 {e e -y] , (A8) 



and that of u e via positron capture on neutrons by 

,2 \ , i A 2/, \ , A 3 



Ql c = ^ n n n e+ [(e e 3 + ) + 3A( e 2 + ) + 3A 2 (e e+ ) + A 3 ] . (A9) 



Here, the number densities of electrons and positrons are described by 



n ^ = 77^3 (^r) 3 ^"2(±%-), (A10) 

(ncy 
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where h is the Planck constant, and the n-th energy moments of electrons and positrons, (e") and (e")*, are given by 



and 



(e:r=(k B T) 



„T2 + n(Ve-^/k B T) 



-F 2 (r?e) 
Pair processes 

The energy emission rates of v e and v e per unit volume via electron-positron pair annihilations are given by 



(All) 
(A12) 



C 2 A-C 2 



36 



8tt 



CqC q 



(A13) 



s (he) 3 J (m e c 2 ) 2 

where Ca = \, CV = 5 + 2 sin 2 # w , and sin 2 6>w = 0.23. Here fpair,^ and are factors for the phase space blocking of v & and v e 

for the pair processes, respectively and the factor for f, is approximately expressed as 

-1 



p . 

1 pair,!/; 



1 +exp 



1 T\(r\ e ) 1 T\(-r) e ) 

+ - — r - Vvi 



2^(r, e ) 2F 3 {-r, e ) 

The energy emission rate of v x per unit volume via the pair annihilations is given by 



(C y -C A ) 2 + (C v + C A -2) 2 / 8tt 
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(hey I (m e c 2 ) 



tor) 9 [J-4(77 C )^3(-77,) + ^3(^)^4(-?7 e )]( J Ppan-,, l ) 2 . 



2^2' 



(A14) 



(A15) 



Nucleon-nucleon bremsstrahlung processes 
The energy emission rate per unit volume through nucleon-nucleon bremstrahlung of a single neutrino pair is well approximated 



with 



28, 



Q% = 1.04 x 10 30 C(Z„ 2 +^ + — X n X p ) 



10 14 gcm" 



5.5 



lMeV 



erg cm 3 s 1 , 



(A16) 



where £ is set to be 0.5 dBurrows et al. 2000), and X„ and X p are the mass fraction of neutrons and protons, respectively. 

Plasmon decay processes 
The emission rates of v e and v e per unit volume via plasmon decay are given by 

cr C 



npl _ dP 1 _ pi ^ 

* D < v 2>a,(hcf (m e c 2 ) 2 ' 



(k B T) S 7 6 e^(l + ~f)P pla , Mc P p l^ r 



(A17) 



where a* = 1/137.036 is the fine-structure constant, 7 = 70 y/r] 2 + 7r 2 /3 with 70 = 2y / a7737r = 5.565 x 10~ 2 , -P P ias,^ and P v \ as .u e 
are factors for the phase space blocking of v e and v e for plasmon decay processes, respectively. The factor for v t is approximately 
expressed as 



plas.fj 



1 +exp 



1 7 2 
1 + -- 



2 1+7 

The emission rate of v x per unit volume via plasmon decay is given by 

47T 3 (TqC 



P?}=(C v -\f 



(fer) 8 7VT(l+7)(Ppi as ,,J Z 



v ' 3a*(/ic) 6 (m e c 2 ) 2 

The energy emission rates of Vj per unit volume via plasmon decaye is given by 



(Ai8) 



(A19) 



(A20) 



Neutrino absorption on helium and inelastic neutrino-helium scatterings 

In addition to the heating processes thro ugh the neutri no absorption on nucleons, the heating processes due to neutrino-helium 
interactions are taken into account, as in Ohnis hTet al.l (120071) . Through the absorption of v e and D e and the neutrino-helium 
inelastic scatterings on nuclei via neutral currents, v + (A,Z) — > v + (A,Z)*, the heating rate per unit volume, Q a , is evaluated as 



Q a = 



pX A 31.6MeV 
m B (r/10 7 cm) 2 



10 52 ergs 



5 Me V \ A" 1 (q+£,e + <^ e A x) r„ e 



T Ve J 



10- 40 cm 2 MeV 



10 52 er ; 



v. ( 5MeV \ A" 1 (<tze Vc + 0% £ e A x ) r± 

rgss" 1 V T% J 10- 40 cm 2 MeV 



10 52 ergs s 



10MeV \ A"' «g£ + <g£) r „ 
T Vx ) 10- 40 cm 2 MeV 



(A21) 
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where is the mass fraction of the nucleus and tob is the atomic mass unit dHaxtonl 19881) . The last term denotes the sum of the 
contributions from fi and r neutrinos. The cross section for each neutral-current is evaluated by the following fitting formula, 



A- 1 (v u v E£ + o u D E& Tu =a 



T„-T 
lOMeV 



(A22) 



where a, (3, an d Tq are given in T able I of iHaxtonl (119881) . and are chosen to be a = 1.24 x 10~ 40 MeV cm 2 , j3 = 3.82, and 
Tq = 2.54 MeV (Ohnishi et al. 2007). In the first and second terms on the right hand side of Eq. ( IA2U . the contributions from the 
charged current reactions, <r+ and a~, are also taken into account according to Table II of lHaxtonl (ll988). that is, A" 1 (er+ E v ^) Tv = 

0.30 x 10" 42 MeVcm 2 for T Vc =4MeV and A" 1 ^ E^) To , = 1.20 x 10" 42 MeVcm 2 for =5MeV. 

The rate of change in specific energy 
The rate of change in energy in unit volume, Qe, is evalutated as 

Qe = (Ql+Q%)+Q a -(Ql+Q%)-lQ b :,-(Qi:+ Qg+4QE)-(Q£+o£+4G£), ( A2 3) 

where a factor of 3 in 3 corresponds three types of a neutrino pair vv. 

THE RATE OF CHANGE IN ELECTRON FRACTION, g N 
The rate of change in Y e per unit volume, 2n, is evalutated as 

g N = (K. - K e -K+ K. ) — , (B i) 

P 

where TZ^ and 1Z% are the absorption rates of i/ e and v e , respectively, and 1Z% and 1Z% are the emission rates of v e and P e through 
the capture of electrons on protons and that of positron on neutrons, respectively. We ignore the variations o f the electron fractio n 
by the neutrino absorption on 4 He, since they are minor and give no qualitative difference to the dynamics (lOhnishi et alJ [2007). 
The absorption rates of v e and D e per baryon are computed as 

L v n n (el )+2A(e„ ) + A 2 , , 
4nr z cf Uc (e„ e ) 

and 

U ^ = ° C 4^cf- <^> ' (B3) 

respectively. 

The emission rate of v e per baryon through the electron capture on protons is given by 

TZl c = l -acn v n s - [< e 2 -)* + 2A(e e ->* + A 2 (e»-)*] (B4) 
and that of P e via the positron capture on neutrons by 

K% = ia C n n « e+ [<ee+)+2A<e e+ ) + A 2 ]. (B5) 



